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Abstract 


The  purpose  of  this  researeh  effort  is  to  experimentally  generate  an  eleetromagnetie 
Gaussian  Sehell-model  beam  from  two  eoherent  linearly  polarized  plane  waves.  The 
approaeh  uses  a  sequenee  of  mutually  eorrelated  random  phase  screens  on  phase-only 
liquid  crystal  spatial  light  modulators  at  the  source  plane.  The  phase  screens  are 
generated  using  a  published  relationship  between  the  screen  parameters  and  the  desired 
electromagnetic  Gaussian  Schell-model  source  parameters.  The  approach  is  verified  by 
comparing  the  experimental  results  with  published  theory  and  numerical  simulation  results. 
This  work  enables  the  design  of  an  electromagnetic  Gaussian  Schell-model  source  with 
prescribed  coherence  and  polarization  properties. 
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EXPERIMENTAL  METHOD  OE  GENERATING  ELECTROMAGNETIC  GAUSSIAN 


SCHELL-MODEL  BEAMS 

I.  Introduction 

Electromagnetic  Gaussian  Schell-model  (EGSM)  beams  use  vector  theory  and  have 
been  proposed  only  recently  within  the  literature  [15].  With  that  said,  EGSM  beams  are 
a  natural  extension  from  Gaussian  Schell-model  (GSM)  beams  which  use  scalar  theory 
[7,  9].  When  considering  EGSM  beams,  interesting  coherence  and  polarization  properties 
have  been  revealed  [21].  During  propagation  a  reduction  in  scintillation  has  been  observed 
as  well  as  changes  in  the  state  of  polarization.  Due  to  these  observations,  EGSM  beams 
have  attracted  special  attention  for  the  potential  use  in  free-space  optical  communications, 
imaging  through  turbulence,  and  remote  sensing  applications  [11,  17,  18].  Improved 
performance  of  the  aforementioned  applications  drives  research  on  this  subject  and  requires 
the  ability  to  control  the  attributes  of  EGSM  sources  [2]. 

Theory  has  greatly  improved  in  understanding  EGSM  beams  over  the  last  few  decades. 
Their  propagation  aspects,  correlation  and  polarization  properties,  and  realizability 
conditions  are  well  studied  and  documented  [8,  16,  24,  25].  Current  research  efforts 
include  analytical  and  experimental  methods  proposed  to  produce  EGSM  sources  [1,  12- 
15,  21,  23,  26].  Many  of  the  proposed  methods  of  generation  are  based  on  interferometer 
designs  and  may  use  either  rotating  phase  wheels,  ground  glass  diffusers,  or  a  more 
practical  design  based  around  a  liquid  crystal  (EC)  spatial  light  modulator  (SLM)  [15,  21]. 
There  are  proposed  methods  and  experimental  designs  for  measuring  the  EGSM  beam 
parameters  [12,  23].  The  most  recent  successful  experiment  demonstrated  a  reduction 
in  scintillation  of  a  completely  unpolarized  EGSM  beam  propagated  through  simulated 
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atmospheric  turbulence  [1].  Practical  techniques  are  all  shown  by  the  above  research 
efforts  to  generate  an  EGSM  beam  and  work  together  to  validate  the  existing  theory; 
however,  these  efforts  do  not  focus  on  the  ability  to  generate  an  EGSM  beam  with  desired 
characteristics. 

1.1  Research  Objective 

The  objective  of  this  research  effort  is  to  design  an  experiment  with  the  capability 
to  generate  an  EGSM  beam  with  prescribed  coherence  and  polarization  properties  from 
two  coherent  linearly  polarized  plane  waves.  The  experimental  design  will  utilize  a  pair 
of  phase-only  nematic  EC  SEMs  which  will  display  a  sequence  of  mutually  correlated 
random  phase  screens.  As  such,  the  desired  EGSM  source  parameters  dictate  the  required 
screen  parameters  [2].  The  random  phase  screens  will  then  be  generated  with  the 
required  parameters  from  the  relationship,  and  the  experimental  results  will  be  verified 
by  comparison  to  published  theory  and  computational  simulation  results  [2,  17]. 

1.2  Limitations 

Many  variables  are  involved  in  the  ability  to  generate  an  EGSM  source  experimentally. 
Most  of  the  limitations  in  this  effort  derive  from  time  allotted  to  complete  the  experiment; 
however,  there  are  limitations  due  to  resources.  Many  of  the  optical  laboratory  resources 
used  in  the  construction  of  the  system  will  be  on-hand  equipment  due  to  the  often 
prohibitive  cost  or  time  necessary  to  acquire  new  equipment.  This  brings  about  limitations 
due  to  age  or  quality  of  the  equipment.  Ultimately  without  more  time,  there  will  be 
equipment  not  calibrated  and  the  full  capability  of  the  experiment  design  will  not  be 
verified.  This  will  be  evident  when  compared  to  results  seen  in  theory  and  simulation. 

1.3  Implications 

Results  from  this  research  effort  could  feasibly  be  applied  directly  to  free  space 
optical  communications,  imaging  through  turbulence,  remote  sensing,  or  directed  energy 
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programs.  If  improvements  in  the  electro-optical  components  and  optical  elements  are 
made,  or  advancements  in  the  use  of  the  SLMs  are  undertaken,  the  experimental  results 
could  be  improved  greatly  allowing  for  further  research  and  possible  application  in  the 
aforementioned  programs. 

1.4  Preview 

The  experimental  research  presented  here  aims  to  demonstrate  that  an  EGSM  source 
can  be  generated  with  the  desired  coherence  and  polarization  properties  with  high 
fidelity  and  configurability.  Chapter  2,  Literature  Review,  discusses  the  knowledge  base 
necessary  for  the  understanding  of  EGSM  sources,  generation,  and  propagation.  Chapter 
3,  Methodology,  details  the  computational  methodology  that  led  to  the  experimental 
design,  as  well  as  the  approach  to  constructing  the  experiment.  Chapter  4,  Analysis  and 
Results,  interprets  the  gathered  theoretical,  simulated,  and  experimental  data.  Chapter 
5,  Conclusions  and  Recommendations,  discusses  the  validity  and  performance  of  the 
experiment  based  on  comparisons  of  the  theoretical,  simulated,  and  experimental  results. 
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II.  Background 


All  optical  fields  undergo  random  fluctuations.  They  may  be  small,  as  in  the  output 
of  many  lasers,  or  they  may  be  appreciably  larger  as  in  light  generated  by  thermal 
sources.  The  underlying  theory  of  fluctuating  optical  fields  is  known  as  coherence  theory. 
An  important  manifestation  of  the  fluctuations  is  the  phenomenon  of  partial  coherence 
and  partial  polarization.  Unlike  usual  treatments  it  describes  optical  fields  in  terms  of 
observable  quantities  and  elucidates  how  such  quantities  change  as  light  propagates  [24]. 
This  background  information  provides  a  unified  treatment  of  the  phenomena  of  coherence 
and  polarization.  Chapter  2  discusses  the  theory  behind  coherence  and  polarization,  which 
are  of  considerable  importance  with  the  propagation  of  EGSM  beams,  as  well  as  discussing 
the  theory  behind  EGSM  beams. 

2.1  Coherence 

In  a  given  random  field,  coherence  describes  the  degree  to  which  one  point  relates  to 
any  other  point  in  the  field  in  time  or  space.  Coherence  is  realized  mathematically  through 
the  correlation  function  r(ri,  r2;  C,  C)-  The  correlation  function  E  depends  on  two  points  in 
space  (ri  and  r2)  or  two  instances  in  time  (C  and  12). 

2.1.1  Mutual  Coherence  Function. 

In  second-order  coherence  theory.  Wolf,  Goodman,  and  others  present  the  mutual 
coherence  function  (MCE)  r(ri,r2,T)  which  is  valuable  in  analyzing  spatial  coherence 
[6,  25].  Specifically, 

r(ri,r2,r)  =  {u{ri,t  +  T)u{r2,t)),  (2.1) 

which  is  a  time  auto-correlation  of  an  analytic  function  u{r,t)  at  two  points  in  space 
(ri  and  r2).  Equation  (2.1)  reduces  to  a  self-coherence  function  when  only  a  single  point  in 
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space  r  is  analyzed.  An  assumption  must  be  made  when  using  Eq.  (2.1)  that  the  field  is  at 
least  wide-sense  stationary,  i.e.,  the  average  field  has  no  explicit  time  dependence  [25]. 
2.1.2  Complex  Degree  of  Coherence. 

The  complex  degree  of  coherence  (CDoC)  y(ri,r2,T)  is  acquired  by  normalizing  the 
MCE  as  shown  in  Eq.  (2.1),  where 


y(ri,r2,T) 


r(ri,r2,T) 

Vr(ri,ri,T)r(r2,r2,r) 


(2.2) 


A  normalized  unit  of  measure  is  given  by  the  magnitude  of  the  CDoC  for  the  amount 
of  temporal  or  spatial  coherence  of  a  field  for  two  points  in  space,  ri  and  r2,  and  a  time 
difference  r  =  t2  -  h.  The  field  is  considered  fully  coherent  if  |y(ri,r2,r)|  =  1,  i.e.,  two 
different  points  in  space  are  correlated.  Conversely,  the  field  is  considered  fully  incoherent 
if  |y(ri,r2,T)|  =  0,  i.e.,  two  points  in  space  are  uncorrelated.  A  field  that  measures 
0  <  |y(ri,r2,T)|  <  1  is  considered  partially  coherent. 

2.1.3  Cross-Spectral  Density. 

The  cross-spectral  density  (CSD)  W(ri,r2,a))  is  another  way  of  analyzing  spatial 

coherence  [25].  The  CSD  and  MCE  form  a  Eourier  transform  pair  given  by 

00 

r(ri,r2,T)  =  ^  J  lT(ri,r2,m)exp(jmT)dT.  (2.3) 

0 

This  demonstrates  the  ability  to  analyze  spatial  coherence  in  the  space-frequency  domain 
with  the  CSD  as  opposed  to  the  space-time  domain  with  the  MCE.  The  following  result, 
derived  by  Wolf  [25], 

Wir  ur2,  CO)  =  {U(ruco)U*(r2,co)),  (2.4) 


shows  the  CSD  is  the  auto-correlation  function  of  an  ensemble  of  sample  functions 
(t/(r,  m)}.  With  the  use  of  the  Wiener- Khinchin  theorem,  the  autocorrelation  and  spectral 
density  form  the  Eourier  transform  pair  for  a  zero-mean,  wide-sense  stationary  random 
process  [6,  25]. 
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2.1.4  Spectral  Degree  of  Coherence. 

Normalizing  the  CSD  as  given  in  Eq.  (2.4),  generates  the  spectral  degree  of  coherence 
(SDoC)  /i(ri,  r2,  oj),  where 


ix{ri,r2,cx)) 


W{ri,r2,<jS) 

^|W(Yl,Yl,^0)W{Y2,r2,C0) 


(2.5) 


A  normalized  unit  of  measure  is  given  by  the  magnitude  of  the  SDoC  for  the  amount  of 
spatial  coherence  of  a  field  for  two  points  in  space,  ri  and  r2,  and  angular  frequency  a». 
Two  different  points  in  space  are  correlated  if  |/i(ri,  r2,  a>)\  =  1  and  the  field  is  considered 
spatially  coherent;  conversely,  two  points  in  space  are  uncorrelated  if  |yu(ri,  r2,  a»)|  =  0  and 
the  field  is  considered  spatially  incoherent.  A  field  that  measures  0  <  |/i(ri,r2,m)|  <  1  is 
considered  spatially  partially  coherent. 

2.1.5  Gaussian  Schell-Model  Source. 

In  the  source  plane  of  the  field  at  the  origin,  the  CSD  is  structured  as 


W{puP2,iO)  =  {U{puiO)U*{p2,CU)), 


(2.6) 


where  pi^  =  xi^2X+yi,2y-  Accordingly,  the  CSD  W(pi,p2,  co)  of  a  GSM  source  is  structured 
as  [25] 

W(pup2,co)  =  ^jS{pi,<j))^jS{p2,oj)p{P2-  Puoj) 

5(p,m)=A2exp|-^j  •  (2.7) 

pip,  co)  =  exp 

Parameters  A^,  w,  and  £  are  space  independent  but  are  dependent  on  angular  frequency  co. 
When  Eq.  (2.7)  is  substituted  into  Eq.  (2.5),  the  magnitude  of  the  SDoC  becomes 

|//(P2 -pi,m)|  =  exp|-i^^^^^^j,  (2.8) 

which  is  only  dependent  on  the  distance  between  two  points  and  not  the  points  themselves 
[25]. 
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The  source  coherence  length  i  is  the  distance  between  two  points  |p2  -  Pi  I  where  the 
magnitude  of  the  SDoC  falls  to  1  of  its  original  on-axis  value.  This  is  a  direct  result 
of  the  relationship  found  in  Eq.  (2.8).  The  GSM  source  is  spatially  coherent  and  the  two 
points  are  correlated  if  |p2  -  Pi  I  ^  however,  the  GSM  source  is  spatially  incoherent  and 
the  two  points  are  uncorrelated  if  |p2  -  pi|  »  i.  The  GSM  source  is  partially  spatially 
coherent  if  0  <  |p2  -  Pil  <  ^. 

2.1.6  Cross-Spectral  Density  Matrix. 

The  cross-spectral  density  matrix  (CSDM)  W(ri,r2,m)  is  utilized  for  the  analysis  of 
spatial  coherence  of  electromagnetic  vector  fields  in  the  space-frequency  domain  [25].  The 
CSDM  is  the  outer  product  generated  from  electric  field  vectors  of  the  following  form: 


E{p,  oj)  =  E^(p,  to)x  -t  Ey(p,  oj)y 

/  N 

EAp,a)) 

Eyip,  CO)^ 


(2.9) 


such  that 

W(pi,p2,m)  ={E(puCo)E^{p2,<^)) 


and 


/  \ 
I  E^ipuco) 

'  [Ey(p2,  co)^ 


(  \ 


E*y{p2,CO)^ 


\e,,ei,)  {e^,e;,) 

SEy.Ei,)  {e^,e;,}^ 


Waf}(Pup2,io)  =  {EJpi,a))E*(p2,a))) 


(2.10) 


(2.11) 


where  ^  denotes  Hermitian  conjugate.  In  Eq.  (2.9),  E^ip,  co)  and  Ey(p,  oj)  are  analytic 
functions  in  two  mutually  orthogonal  directions  perpendicular  to  the  direction  of 
propagation. 
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Accordingly,  the  SDoC  is  determined  from  the  CSDM  using  the  following  relationship 


[25]: 


,  .  Tr{W(pi,p2,t^)} 

iu(pi,p2,a))  =  ,  . — 

^Tr  { W(pi ,puco)}  ^JTr:  { W(p2,  P2,  oj)] 


(2.12) 


where  Tr{---}  denotes  the  trace.  Formulated  with  electromagnetic  vector  fields,  a 
normalized  unit  of  measure  is  given  by  the  magnitude  of  the  SDoC  for  the  degree  of  spatial 
coherence,  i.e.,  0  <  \p(pi,p2,(i>)\  <  1. 

The  CSDM  of  a  GSM  source  takes  the  following  element-based  form  [25]: 


Wa/3(fil,P2,iO)  -  sjS  a(Pl,Oj)  ^JS  y3(p2,  0))p{p2  “Pi,  n>) 

IPP 


Sm(p,io)  =  A^exp 


PapiPl  -P\,0J)  =  Baj}  exp 


2wl 


‘  \p2  -p\?^ 

K,  J 


(2.13) 


(a  =  x,y  and  /3  =  x,y) 


Element-based  parameters  A^,  Bap,  Wa,  and  tap  are  space  independent  but  are  dependent  on 
angular  frequency  co. 


2.2  Polarization 

Given  electromagnetic  vector  fields  and  the  CSDM  as  defined  in  Eq.  (2.10),  the 
following  relationships  of  interest  exist  for  polarization.  The  first  relationship  of  interest 
is  the  space-  and  angular-frequency-dependent  degree  of  polarization  (DoP)  P{p,  00)  [25] . 
Specifically, 


Pip,  CO) 


4Det{W(p,p,  m)} 


1- 


(2.14) 


(Tr{W(p,p,  m)})^’ 

where,  Det  {•  •  • }  denotes  the  determinant  operation.  A  normalized  unit  of  measure  is  given 
by  the  DoP  for  the  amount  of  polarization  in  a  field  [25].  The  field  is  polarized  when 
P(p,(o)  =  1;  conversely,  the  field  is  unpolarized  when  Pip,co)  =  0.  The  field  is  partially 
polarized  when  0  <  Pip,  co)  <  1. 
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The  last  relationship  of  interest  is  the  single  point  Stokes  vector  Si(p,a)),  where 


I  =  1, 2, 3, 4.  The  components  of  this  vector  are 


S o(p,  m)  =W^Ap, P,  co)  +  Wyyip, P,  oj) 


Si(p,co)  =W^Ap,P,  co)  -  Wyy(p,p,  O)) 


S2(p,co)  =Wj,y(p,p,co)  +  Wy^(p,p,a)) 

S  3(p,  0))  =  j  [WyAp,p,  OJ)  -  W^y(p,p,  m)] 


(2.15) 


The  Stokes  vector  is  useful  in  analyzing  polarization  [25].  Utilizing  the  Stokes  vector,  the 
DoP  is 


P(p,  CO)  = 


\{p,  <jO)  +  S  lip,  <jO)  +  S  lip,  <j0) 


S  o(p,  co) 

which  demonstrates  a  second  method  within  the  analysis. 


(2.16) 


2.3  Electromagnetic  Gaussian  Schell-Model  Beam  Theory 

Now  that  the  fundamental  theory  has  been  covered,  the  theoretical  solution  developed 
by  Korotkova  and  Hyde  for  the  generation  of  the  EGSM  beam  is  reviewed  [2,  17]. 
Korotkova’s  method  uses  the  tensor  technique  to  characterize  the  source.  The  following 
analysis  presents  the  key  equations  from  Korotkova’s  derivation  that  are  necessary  for 
analytical  comparison  of  the  EGSM  source  plane  [17]. 

Beginning  in  the  source  plane,  the  tensor  notation  for  each  of  the  four  CSDM  elements 
of  the  EGSM  have  the  form 


Wapir,  0)  =AaApBapexp 


(a  =  x,y  and  /3  =  x,y) 


(2.17) 


where  f  is  the  4x4  vector  such  that  f  =  (ri,  r2),  ri  and  r2  are  two-dimensional  vectors  in 
the  source  plane,  k  is  the  wave  number,  and  is  the  4x4  matrix  expressed  as 


ik\2<rl  +  ^ 


—  I 

Tx2  -*■ 


—  T  T 

A  1 20-2  Sl^ 


(2.18) 
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Parameters  A^,  Bap,  cr^,,  and  6ap  in  Eq.  (2.17)  are  independent  of  position  but  depend  on 
frequency,  and  I  is  the  identity  matrix. 

The  above  detailed  equations  provide  the  theoretical  solution  for  comparison  of  an 
EGSM  beam  in  the  source  plane.  Now  that  the  source  plane  is  derived,  the  theoretical 
solution  by  Hyde  for  an  EGSM  beam  propagated  through  a  lens  to  the  observation  plane  is 
reviewed.  The  following  are  the  key  equations. 

Each  of  the  four  CSDM  elements  of  the  EGSM  beam  in  the  source  plane  have  the 


form 


Wap{p'i,p'2,^  )  =AaApBapC,X^ 

1 


[  pfl 

^  ip;-P2p^ 

exp 

4(7-2 

exp 

J 

Wap(pi,P2,  f)  -AaApB, 


ap 


exp 


w 

\p\  -  p'l? 


exp 


k,Pi-p^) 


Ifff 


exp 


4cr^ 


exp 


.ft 


exp 


jk 
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(P1P1-P2P2)  dVidV2 


(a  =  x,y  and  =  x,y) 


(2.19) 


where  z  =  0“  is  the  source  plane  directly  before  the  lens  and  z  =  /  is  the  focal  plane  of  the 
lens,  i.e.,  the  observation  plane. 

After  evaluating  the  integrals,  the  four  CSDM  elements  of  the  EGSM  beam  are 


Wap(p\,P2,f) 


AaApBapTT^ 


1 


^a,ap^p,ap  ^ap 

\  (> 

exp 

exp 


-exp 


exp 


'a/S,apxl  -  2bapXiX2  +  aa,apxl\ 

4/2 

k2 

'  <^P.afiy\  -  '^bapyiyi  +  aa,a(iy\\ 

4/2 

(2.20) 


(a  =  x,y  and  =  x,y) 
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where 


and 


4cr^ 


ap,ap  2  +  2^2  ’ 


P?  =(2f?+y?) 
p\  =  (^2+y^) 


(2.21) 


(2.22) 


The  above  detailed  equations  provide  the  theoretical  solution  for  comparison  of  an  EGSM 


beam  in  the  observation  plane. 
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III.  Methodology 


Chapter  3  discusses  a  published  method  of  numerically  generating  an  EGSM  source 
with  prescribed  coherence  and  polarization  properties  and  the  relationships  between  the 
EGSM  source  and  phase  screen  parameters.  These  relationships  provide  the  basis  for 
creating  the  phase  screens  necessary  to  experimentally  generate  an  EGSM  source.  Eurther, 
the  experimental  design  will  be  discussed  in  detail,  to  include  the  calibration  of  the  SEMs 
and  the  collection  and  measurement  of  data. 

3.1  Electromagnetic  Gaussian  Schell-Model  Source  Generation 

With  a  theoretical  solution  for  comparison,  the  method  for  creating  an  EGSM  source 
must  be  detailed.  The  following  method  to  create  the  source  was  developed  by  Hyde  and 
the  key  equations,  steps,  and  much  of  the  derivation  are  reproduced  in  the  following  section 
[2].  The  following  analysis  uses  the  EGSM  source  CSDM  from  Wolf  as  shown  in  Eq.  (2.13) 
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and  reproduced  here  for  convenience,  i.e.,  [25]: 


^(7(Pi’P2,0,  cu)  -  ^JS  iipuCij)  a/s  jip2\ co) 
A’XIPi  - P2I; (^)  O’  =  2c,y  j  =  x,y) 


Si(p-,co)  =A^  exp 

f  -\Pi  -  P2?' 

Piji\p\  -  P2I;  oj)  =  Bij  exp  - - 


Bij  =1  i  =  j 
\Bij\  <1  ii^  j 


Bij  =  B*j 


^ij  =  ^Ji 


1  1 

The  next  step  in  the  analysis  is  to  define  the  electric  field  [reference  Eq.  (2.9)]  in  the 
source  plane,  i.e.,  z  =  0,  as 

E^(p)  =  Cacxp  j  exp  [)0Q,(p)]  (or  =  x,y),  (3.2) 

where  Cq,  =  \Ca\  exp(j  6^)  is  a  complex  constant.  Performing  the  autocorrelations  necessary 
to  fill  the  CSDM  in  Eq.  (2.10)  produces 


2  2 
PI  ,  P2 


(EaiE^i)  =  “I  4)3 +  4^  I  <exp[)<^„(pi)]exp[-j0;j(p2)]> 


a  =  x,y 

^  =  2c,y 


where  cpa  and  (pp  are  random  phase  screens.  These  phase  screens  are  sample  functions 
drawn  from  two  correlated  Gaussian  random  processes  which  are  detailed  in  the  next 


section. 


To  allow  for  the  approximation  of  the  normalized  cross-correlation  function,  the 


function  is  taken  to  be  Gaussian-shaped  and  the  standard  deviations  of  the  phase  screens. 
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and  cr^^,  are  assumed  to  be  greater  than  n.  This  yields 


<exp[i^„i]exp[-j0^2]> 


1  /  2  o  ,  2  \ 

IPl  -P2P 

exp 

“2  V  <l>a  ~  ^P<Pa>Pli^<Pa^>PlS  ^ 

exp 

! P‘l>a‘l’P  ^ h  __ 

where  0  <  <  1  is  a  correlation  coefficient  =  1  if  a  =  yS).  Eq.  (3.4)  is  substituted 

into  Eq.  (3.3)  and  simplified  to  form 


{EaiE*2)  =CQ,C^exp 


P' 


i  P\  ,  Pi  \ 

[  (4^^  4cr2) 

exp 

“9  “  '^P‘i>c‘i>p^‘i>c^‘i>p + 


exp 


IPi  -  P2P 


(3.5) 


^<l>a(pp  / P‘l>a‘l>P  ^ <!> 

The  “self’  terms  of  the  CSDM  are  created  by  letting  a  =  fi  =  x  ox  y,  shown  for  x  as 
follows: 


{E^,e:^)  =  IC.pexp 


/p?+p2y 

IPi  -  P2P 

\  4(r2  ) 

exp 

-  ^<l>x<l>xl^<l>x  - 

(3.6) 


Comparing  Eq.  (3.6)  (or  the  equivalent  for  a  =  yS  =  y)  to  Wolf’s  GSM  form,  the  following 
relationships  are  required: 

1 


5, 

6. 


V2 

1 


|C, 


|C,|=A, 


(3.7) 


Similarly,  the  “cross”  terms  of  the  CSDM  are  produced  by  letting  a  =  x  and  /3  =  y,  from 
which  the  following  relationships  are  required  to  match  Wolf’s  GSM  form 

1 


■'xy 


V2  yP4>x(l>y^<l>x^^y 


Bxy\  =  exp 


■9  [^ix  ~  ^P<l>x^y^<l>x(^>Py  + 


(3.8) 


^Bjcy  =  Ox  -  Oy 

Eurther,  letting  a  =  y  and  p  =  x  provides  the  complement  relationship  to  Eq.  (3.8),  which 
satisfies  Wolf’s  requirement  for  Bij  and  6ij  as  defined  in  Eq.  (3.1).  As  demonstrated  in 
further  detail  in  Ref.  [2],  the  above  approach  for  simulating  an  EGSM  source  is  analytically 
sound. 
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3.2  Phase  Screen  Generation 

Now  that  the  EGSM  source  has  been  created,  the  correlated  random  phase  screens 
must  be  generated.  As  with  the  previous  section  for  creating  a  source,  the  following  method 
for  creating  the  phase  screens  was  completed  by  Hyde  and  the  key  equations  and  steps  are 
highlighted  [2]. 

Let  (p  and  0  be  Fourier  transform  pairs,  i.e.. 


oo 

ff^u. 


y)exp(-j27rfj,x)exp(-j27rfyy)dxdy 


OO 


(3.9) 


(p{x,y)=  (l){f^,fy)&xp{-]2nf^x)&xp{-]2nfyy)dfAfy 


It  must  also  be  noted  that  the  phase  screens  are  zero  mean  and  Gaussian  correlated: 


{(pAx,y))  =  <0>.(v:,y)>  =  {<pa{x,y))  =  0 
{(pa{xi,yi)(pl{x2,y2))  =  (^^exp 


V  ) 


(3.10) 


Expanding  (pa  in  a  Fourier  series  and  taking  the  autocorrelation  produces 


{(pa{Xi,yi)(p*a{X2,y2))  =  EE 


m,n  p,q 


i^imxi  -px2) 


exp 


]^inyi-qy2) 


(3.11) 


where  the  Fourier  series  coefficients  cpamn  and  tpapq  are  zero  mean  Gaussian  random 
numbers  and  L  =  AA  is  the  size  of  the  discrete  grid.  This  expression  must  equal  the 
autocorrelation  of  (pa  computed  using  Eq.  (3.9),  that  is, 

oo 

{(pa{xuy\)(p*a{x2,y2))  =  JJ fy)^^P  [j^^fx(xi  -  X2)]  cxp  [j2;r/j.(yi  -  y2)]  df^dfy, 

—  OO 

(3.12) 

where  is  equivalent  to  the  power  spectral  density  (PSD)  of  (pa'- 


%.^SfxJy)  =  <^4^„exp  +  f^)] . 


(3.13) 
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The  correlation  in  Eq.  (3.12)  must  be  discretized.  To  do  so  the  integrals  must  be 
expanded  in  Riemann  sums.  These  sums  when  compared  to  Eq.  (3.11)  result  in  the 
following  relationships: 


amn^( 

<l‘^. 


apij)  ^<f>a<Pa 
">  = 


m  n\  1 

m  n\  \ 

T'ul? 


(3.14) 


where  {\(pamn?')  is  equivalent  to  the  variance  of  the  Eourier  series  coefficients  (pamn  and  5mp 
and  6nq  are  Kronecker  deltas. 

Thus,  the  phase  screen  (pa  can  be  produced  by  generating  a  matrix  of  unit  variance 
circular  complex  Gaussian  random  numbers  multiplying  Va  by  the  square  root  of 
Eq.  (3.14),  and  performing  a  two-dimensional  discrete  inverse  Eourier  transform,  namely. 


(3.15) 


(pa\i,j]  =  'Yjra[m,n\ 


-exp 


I  m  (  n 


l.ln  \  Lin  \ 
explj— milexplj— njl 


The  output  of  the  inverse  Eourier  transform  Eq.  (3. 15)  is  a  complex  matrix  where  either  the 
real  or  imaginary  part  can  be  used  to  create  (pa.  Here,  the  real  part  is  used. 

To  simulate  the  “cross”  terms  of  the  CSDM,  the  cross-correlation  of  Eq.  (3.15)  must  be 
computed.  Making  use  of  common  trigonometric  identities  and  Euler’s  formula,  additional 
simplifications  can  be  made  to  the  resultant  yielding 


(Re(0;,[/,  j])  Re((py[k,  /]))  =  E  Ip ^  ^(^^Ipxlpx  ^<Py  (py  ) 


exp  <  -n 


In 


(f-  +/2  N 

'  ^(px^x  ^  ^My  ' 


m 

NA 


(3.16) 


N 


exp -  k)  exp \i—n(nj  -  1) 


In 


N 


1 


{NAf 


By  comparing  the  discrete  function  being  transformed  in  Eq.  (3.16)  to  the  continuous 
cross-power  spectral  density  function,  i.e.. 


^.Px^yifxJy)  =  (r^,cr^^,7rp^,^/J^^^exp  [- 


■^^^IxiPy^fx  + 


fy)]  ’ 


(3.17) 
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the  following  relationships  are  obtained: 


P<t>x<l>y 


^  ^  P^xh  .  (3.18) 

'^^'l>x<l>x^lpy<f>y 

_  '^^‘Px4>x^'Py4>y 
P<l>x‘l>y  p2  I  p2 

^:l>x<Px  ^<l>y<l>y 

As  demonstrated  in  further  detail  in  Ref.  [2],  this  method  creates  correlated  random 
Gaussian  phase  screens  necessary  for  the  simulation  and  experimental  generation  of  an 
EGSM  source  to  control  degrees  of  coherence  and  polarization. 


3.3  Spatial  Light  Modulators 

In  this  experiment,  a  dual  SLM  design  serves  as  the  active  controller  of  the  EGSM 
beam  by  means  of  displaying  the  random  phase  screens.  Erom  an  understanding  of  the 
structure  and  operation  of  the  SEM,  potential  sources  of  error  can  be  identified.  This 
section  begins  with  descriptions  of  the  relevant  terminology,  followed  by  discussions  of 
how  inherent  design  and  manufacturing  flaws  are  calibrated  out. 

3.3.1  Principles  of  Operation. 

There  are  several  different  types  of  SEMs  available  for  use,  providing  choices  between 
transmissive  and  reflective  SEMs,  which  may  control  either  phase,  amplitude,  or  both. 
SEMs  have  an  expanding  role  in  several  optical  areas  where  light  control  on  a  pixel- 
by-pixel  basis  is  critical  for  optimum  system  performance.  The  SEMs  chosen  for  this 
experiment  are  Boulder  Nonlinear  Systems  (BNS)  Model  P5 12-0635  XY  Series  EC  SEMs 
with  a  512  X  512  pixel  array  and  15//m  pitch,  thus  the  focus  of  the  following  discussion 
will  be  based  on  electrically  addressed  phase-only  nematic  EC  SEMs  [3]. 

The  structure  of  a  liquid  crystal  on  silicon  (ECoS)  SEM  has  several  defining 
characteristics.  SEMs  control  light  based  on  a  fixed  spatial  (pixel)  pattern.  The  spacing 
between  the  centers  of  pixels  in  the  pattern  is  referred  to  as  pixel  pitch.  By  design,  polarized 
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light  enters  the  device  passing  through  a  cover  glass,  transparent  electrode,  and  LC  layer. 
Beneath  the  LC  layer  are  reflective  aluminum  pixel  electrodes.  The  light  reflects  off  this 
electrode  layer  and  returns  on  the  same  path.  A  voltage  induced  electric  field  between 
the  pixel  electrode  and  the  transparent  electrode  on  the  cover  glass  changes  the  optical 
properties  of  the  LC  layer.  Because  each  pixel  is  independently  controlled,  a  phase  pattern 
may  be  generated  by  loading  different  voltages  onto  each  pixel  [3]. 

The  chosen  BNS  SLMs  are  optimized  to  provide  a  full  wave  (27r  rad)  of  phase  stroke 
upon  reflection  at  the  A  =  635nm  wavelength.  These  SLMs  only  provide  phase  modulation 
when  the  input  light  source  is  linearly  polarized  along  the  vertical  axis.  Additionally,  the 
reflective  pixel  structure  associated  with  a  LCoS  SLM  backplane  acts  as  an  amplitude 
grating  that  diffracts  some  light  into  higher  orders.  In  the  experimental  design,  the  SLMs 
are  aligned  to  utilize  reflected  light  off  of  the  first  diffraction  order  rather  than  the  zeroth- 
order.  An  eight-step  pixel  grating  is  applied  to  the  SLMs  to  direct  energy  into  the  first 
diffraction  order.  The  rationale  for  this  adjustment  will  be  addressed  in  a  later  section. 

3.3.2  Phase  Response  Calibration. 

A  SLM  has  a  unique  response  in  converting  phase  to  digital  command  values. 
Several  methods  for  calibrating  phase  response  were  considered,  such  as  placing  the 
SLM  in  a  Michelson  interferometer,  double-slit  aperture  method,  and  using  amplitude 
modulation  [19].  A  diffractive  amplitude  modulation  method  was  chosen  because,  unlike 
in  interferometry,  the  measurements  are  insensitive  to  vibrations. 

To  perform  the  amplitude  modulation,  the  SLMs  were  aligned  to  use  light  reflected 
into  the  zeroth-order.  The  reflected  light  then  passes  through  a  lens  after  each  SLM  as 
shown  in  Fig.  3.1.  An  imaging  sensor  is  placed  at  the  respective  focal  plane  of  each  lens. 
Prior  to  the  light  entering  the  detector,  it  passes  through  a  quarter-wave  plate  (QWP)  and 
linear  polarizer  (LP),  which  are  required  to  collect  the  necessary  irradiance  to  calculate  the 
phase  response.  This  method  is  thoroughly  documented  by  Schmidt  [19]. 
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The  SLMs  interface  with  a  controller  providing  8-bits  of  pixel  data.  An  array  of 
commands  is  loaded  into  the  SLM  memory.  Each  command  from  0  to  255  is  a  flat  phase 
image,  i.e.,  every  addressable  pixel  in  the  array  is  set  to  the  same  value  for  each  command. 
The  LP  prior  to  each  SLM  is  aligned  to  45-degrees,  the  QWP  prior  to  the  camera  is  aligned 
to  0-degrees,  and  the  LP  prior  to  the  camera  is  aligned  to  45-degrees  for  the  first  data 
collection.  Each  command  is  stepped  through  on  each  SLM  and  the  irradiance  on  the 
sensor  is  measured.  Using  only  this  first  measurement  allows  the  phase  to  be  computed 
over  only  half  of  the  unit  circle  [19].  To  complete  the  second  measurement,  the  LP  prior 
to  the  SLM  remains  aligned  to  45-degrees,  the  QWP  is  then  aligned  to  45-degrees,  and  the 
LP  remains  aligned  to  45 -degrees.  Each  command  is  again  stepped  through  on  each  SLM 
and  the  irradiance  is  measured.  With  the  second  measurement,  phase  can  be  measured  over 
the  entire  unit  circle,  and  a  standard  unwrapping  technique  can  be  used  to  compute  the  true 
physical  phase  commanded  to  the  SLM  [19]. 


SLM 
Path  2 


Lens  QWP  Camera 


Ligure  3.1:  Schematic  illustration  of  system  setup  for  linear  phase  response  calibration. 
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The  irradiance  measurements  for  eaeh  SLM  with  and  without  the  QWP  are  plotted 
generating  a  power  eurve  with  maximums  and  minimums  where  eaeh  phase  step  is  equal  to 
n  radians  of  phase,  similar  to  Fig.  3.2(a)  and  Fig.  3.3(a).  To  obtain  the  phase,  the  irradianee 
data  is  input  into  a  four-quadrant  arctangent  funetion 


^^SLM  -  tan 


/,  -r 


(3.19) 


where  h  and  h  are  the  irradianee  data  eolleeted  with  the  two  measurements.  This  ealeulates 
phase  that  is  wrapped  over  the  interval  {-n,  :7r].  This  wrapped  phase  ean  then  be  unwrapped 
[19].  hln  region  ean  then  be  seleeted  to  generate  a  look-up-table  (LUT)  whieh  maps  the 
desired  phase  eommand  to  the  eleetrieal  eommand  for  eaeh  pixel. 

During  this  ealibration,  ditheulty  was  experieneed  generating  power  eurves  without 
saturation.  The  manufaeturer  provides  linear  phase  response  LUTs  for  eaeh  SLM.  To  verify 
the  SLMs  were  working  properly,  the  manufaeturer  LUTs  were  loaded  to  eaeh  SLM  and 
the  same  proeedure  above  was  performed,  stepping  through  eaeh  eommand.  What  would 
be  expeeted  if  the  LUTs  were  an  aeeurate  ealibration  would  be  a  near  linear  unwrapped 
eommand  to  phase  plot  of  a  2n  region.  As  shown  in  Fig.  3.2(b)  and  Fig.  3.3(b),  the 
manufaeturer  provided  LUTs  were  adequate  and  used  throughout  this  researeh  effort. 
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SLM  1  SLM  1 


Figure  3.2:  Linear  phase  response  calibration  plots  for  SLM  1  using  manufacturer  provided 
look-up-table  showing  (a)  measured  irradiance  with  and  without  a  quarter-wave  plate  and 
(b)  the  unwrapped  phase. 
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SLM2 


SLM2 


Figure  3.3:  Linear  phase  response  calibration  plots  for  SLM  2  using  manufacturer  provided 
look-up-table  showing  (a)  measured  irradiance  with  and  without  a  quarter-wave  plate  and 
(b)  the  unwrapped  phase. 


To  smooth  the  curves  and  optimize  the  speed  of  converting  phase  to  commands,  a  Fade 
fit  is  applied  to  the  LUT  data  for  each  SLM  as  shown  in  Fig.  3.4.  A  Fade  function,  P{x), 
has  the  form 

M 

P{x)  =  - ,  (3.20) 

1  +  Z 

where  where  F  is  the  dependent  variable,  x  is  the  independent  variable,  m  and  n  are  integer 
indices,  M  and  N  are  the  highest  polynomial  orders  in  the  numerator  and  denominator, 
respectively,  and  where  and  bn  are  the  polynomial  coefficients  in  the  numerator  and 
denominator,  respectively. 
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To  use  this  calibration  for  commanding  phases  onto  the  SLM,  phase  values  (one  value 
for  each  pixel)  on  the  interval  [0,  Itt)  are  sent  to  the  SLM  and  converted  directly  to  an 
array  of  command  values  using  a  Fade  function  with  the  appropriate  coefficients  rather 
than  searching  a  LUT  [19].  This  greatly  improves  the  performance  and  efficiency  of  the 
SLM  control  code. 


SLM  1 


Figure  3.4:  Fade  fit  to  manufacturer  provided  look-up-table  data  for  SLM  1  and  SLM  2. 


3.3.3  Static  Aberration  Calibration. 

A  LCoS  SLM  has  an  inherent  static  aberration  across  its  full  aperture  due  to  the 
manufacturing  process  [19].  Thus,  the  shape  of  this  aberration  must  be  measured  to  apply 
the  appropriate  phase  map  to  the  SLM  to  compensate.  There  are  many  different  methods 
to  perform  a  static  aberration  calibration. 
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Parametric  optimization  was  the  first  method  considered  but  was  prohibited  by  the 
requirement  to  capture  an  image  in  both  the  lens  plane  and  focal  plane.  The  SLM  aperture 
is  too  large  to  fully  fit  on  the  sensor  array  and  it  was  not  possible  to  demagnify  the  SLM 
aperture  given  the  testbed  and  limited  amount  of  space.  The  next  calibration  method 
considered  used  interferometry,  which  provides  a  direct  measurement  of  the  SLMs  inherent 
phase  distortion.  The  measurement  comes  from  interferogram  analysis  using  a  Michelson 
interferometer  setup.  This  method  was  attempted  with  little  success.  An  interferogram  was 
somewhat  visible  but  this  method  was  ultimately  too  sensitive  to  vibrations. 

The  last  method  considered  and  chosen  was  an  iterative  Fourier  transform  method, 
using  forward  and  inverse  Fourier  transforms  to  propagate  a  field  back  and  fourth  between 
the  pupil  and  focal  planes  while  imposing  physical  constraints  at  both  planes.  This 
approach  to  phase  retrieval  is  known  as  the  Gerchberg-Saxton  algorithm  [10]. 

The  Gerchberg-Saxton  algorithm  is  an  error-reduction  algorithm,  as  the  error  in  the 
solution  decreases  with  every  iteration  [10].  Convergence  occurs  when  the  decrease  in  error 
stagnates.  Phase  diversity  is  used  here  to  extend  the  original  Gerchberg-Saxton  algorithm 
[10].  Figure  3.5  is  a  flow  chart  showing  the  modified  algorithm  with  phase  diversity. 
Stepping  through  the  chart,  an  initial  guess  for  the  static  aberration  is  made.  The  phase 
function  is  multiplied  by  the  pupil  function.  Then,  a  known  phase  aberration  is  added  to 
the  static  phase  estimate.  Next,  a  Fourier  transform  translates  the  pupil  field  to  the  focal 
plane.  Here  the  calculated  amplitude  is  substituted  by  the  measured  amplitude  from  the 
known  phase  aberration.  An  inverse  Fourier  transform  then  translates  the  modified  field 
to  the  pupil  plane  where  the  applied  known  aberration  is  subtracted  from  the  resulting 
phase.  This  constitutes  one  iteration;  the  remaining  phase  is  an  estimate  for  both  the 
pupil  amplitude  and  static  aberration  [10].  Each  successive  iteration  through  this  loop 
ideally  uses  a  new  pair  of  known  pupil  phase  and  measured  focal  plane  amplitudes  with  the 
static  phase  aberration  estimates  improving  each  time.  For  practicality,  the  known  phase 
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aberrations  and  captured  amplitudes  are  a  finite  set  of  measurements  and  the  algorithm 
loops  through  them  until  a  stopping  criterion  is  met.  What  remains  when  the  algorithm 
stops  is  the  static  aberration  phase  estimate. 


Figure  3.5:  Modified  Gerchberg-Saxton  algorithm  flow  chart  for  iterative  Fourier  transform 
calibration  method. 


After  successfully  verifying  the  algorithm  in  simulation  using  Zernike  polynomials, 
the  algorithm  was  ready  to  be  used  to  calibrate  the  SLMs.  The  experimental  procedure 
first  required  the  generation  of  unique  phase  screens  to  be  applied  to  each  SLM.  To 
provide  phase  diversity,  Zernike  terms  3-9  were  used  to  make  three  sets  of  phase  screens 
with  varying  peak  aberration  coefficients,  totaling  21  phase  screens.  These  screens  were 
commanded  to  each  SLM  and  the  corresponding  diffraction  patterns  were  gathered  at 
the  focal  plane.  The  known  phase  screens  and  corresponding  collected  images  were 
then  loaded  into  the  iterative  Fourier  transform  algorithm  to  loop  through  the  images 
as  described  above  with  a  set  maximum  of  300  iterations  or  until  a  sum-squared-error 
threshold  was  reached.  Each  SLM  was  calibrated  individually  as  shown  in  Figs.  3.6  and 
3.7. 
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Figure  3.6:  Schematic  illustration  of  system  setup  for  static  phase  aberration  calibration  for 
Path  1. 


Iris 


Figure  3.7:  Schematic  illustration  of  system  setup  for  static  phase  aberration  calibration  for 
Path  2. 
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3.3.4  Comparison  of  Results. 

The  intention  of  performing  both  the  linear  phase  response  and  the  static  aberration 
calibrations  is  to  flatten  the  wavefront  at  the  pupil  plane,  leaving  only  the  effect  of 
diffraction  by  the  pupil  itself.  The  square  aperture  of  the  SLM  is  the  pupil,  the  far-field 
diffraction  pattern  is  expected  to  be  a  2-D  sine  pattern.  The  static  aberration  estimates 
gathered  from  the  above  method  for  each  SLM  are  shown  in  Fig.  3.8.  The  observed 
diffraction  pattern  for  each  SLM  are  shown  without  correction  in  Figs.  3.9(a)  and  3.10(a) 
and  with  correction  in  Figs.  3.9(b)  and  3.10(b). 

When  performing  the  static  aberration  calibration,  the  initial  alignment  of  the  SLMs 
used  the  light  reflected  into  the  zeroth-order  as  described  in  the  linear  phase  response 
calibration.  The  static  calibration  was  repeatedly  failing  to  provide  a  good  aberration 
estimate.  It  was  eventually  discovered  there  was  a  specular  reflection  off  the  front  face 
of  SLM  glass  that  was  not  being  controlled  by  the  SLMs.  This  reflection  traveled  the  same 
optical  path  to  the  sensor  and  was  placing  a  bright  spot  in  the  amplitude  measurements  for 
which  the  algorithm  could  not  account.  To  correct  this  issue,  the  alignment  of  the  SLMs 
was  adjusted  to  use  the  light  reflected  into  the  first-order  as  stated  above.  This  adjustment 
was  not  made  for  the  linear  phase  response  calibration  and  it  is  unknown  if  this  anomaly 
affected  those  calibration  results. 
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SLM  1  Static  Aberration  Estimate 
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SLM  2  Static  Aberration  Estimate 
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Figure  3.8:  Static  phase  aberration  estimates  for  (a)  SLM  1  and  (b)  SLM  2. 
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SLM  1  without  Aberration  Correction 
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Figure  3.9:  Static  aberration  calibration  diffraction  patterns  showing  a  flat  phase  applied  to 

SLM  1  and  propagated  through  the  system  to  the  observation  plane  (a)  without  correction 

applied  for  the  static  aberration  and  (b)  with  correction  applied  for  the  static  aberration. 
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SLM  2  without  Aberration  Correction 
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SLM  2  with  Aberration  Correction 
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Figure  3.10:  Static  aberration  calibration  diffraction  patterns  showing  a  flat  phase  applied  to 

SLM  2  and  propagated  through  the  system  to  the  observation  plane  (a)  without  correction 

applied  for  the  static  aberration  and  (b)  with  correction  applied  for  the  static  aberration. 
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3.4  System  Model 

3.4.1  Design. 

For  the  experiment  designed  for  this  researeh,  the  source  chosen  to  generate  the  EGSM 
beam  is  a  HeNe  gas  laser,  whose  radiation  is  almost  completely  coherent  and  completely 
polarized,  with  an  output  wavelength  of  d  =  632. 8nm.  As  shown  in  Fig.  3.11,  the  beam 
from  the  laser  source  passes  through  a  beam  expander.  After  leaving  the  beam  expander,  the 
expanded  beam  passes  through  an  iris.  The  beam  expander  is  adjusted  to  fill  a  minimum 
region  of  interest  while  also  maintaining  collimation.  This  was  tested  after  each  optical 
element  discussed  below.  The  iris  was  adjusted  to  prevent  the  beam  from  over-filling  the 
optical  elements  while  still  filling  the  region  of  interest. 
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Figure  3.11:  Schematic  illustration  of  experiment  design  to  generate  an  EGSM  source. 
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After  the  iris,  the  beam  passes  through  a  half-wave  plate  (HWP)  and  LP.  The 
purpose  of  these  optical  elements  is  to  control  the  amplitude  of  the  beam.  Following 
these  elements,  the  beam  enters  a  polarizing  beam  splitter  (PBS).  The  PBS  allows  the 
horizontally  polarized  light  to  pass  through  and  reflects  the  vertically  polarized  light.  These 
beam  components  constitute  the  Path  1  (horizontally  polarized  axis)  and  Path  2  (vertically 
polarized  axis),  as  referenced  from  this  point  forward  and  labeled  in  Fig.  3.1 1. 

As  previously  discussed  in  Section  3.4,  the  SLMs  used  in  this  experiment  act  only  on 
vertically  polarized  light.  The  portion  of  the  beam  passing  through  the  PBS  is  horizontally 
polarized.  After  leaving  the  PBS,  this  component  passes  through  another  HWP  and  LP.  The 
HWP  is  used  to  control  the  relative  amplitude  {Ay)  of  Path  2  of  the  EGSM  beam.  The  LP 
is  used  to  ensure  only  vertically  polarized  light  is  incident  on  the  SLM.  The  portion  of  the 
beam  reflected  by  the  PBS  is  vertically  polarized.  After  leaving  the  PBS,  this  component 
also  passes  through  another  HWP  and  LP.  Again,  the  HWP  is  used  to  control  the  relative 
amplitude  {A^}  of  Path  1  of  the  EGSM  beam.  The  EP  is  used  to  ensure  only  vertically 
polarized  light  is  incident  on  the  SEM. 

Both  beam  components  are  now  incident  on  the  EC  SEMs.  The  SEMs  display  the 
correlated  Gaussian  random  phase  screens.  The  SEMs  are  placed  and  oriented  to  reflect 
the  beams  parallel  to  each  other,  with  a  grating  applied  to  direct  the  beam  energy  into  the 
first-order. 

Eollowing  the  SEM  on  Path  7  is  a  4-/  lens  system  (ESI)  shown  in  Eig.  3.12  which 
serves  multiple  purposes.  This  lens  system  has  an  iris  at  the  focal  plane  between  the 
lenses.  Both  lenses  used  are  plano-convex  and  have  a  350mm  focal  length  so  as  to  not 
magnify  or  demagnify  the  beam.  The  first  purpose  of  the  lens  system  is  to  remove  unwanted 
diffraction  orders  with  the  iris.  If  the  diffraction  orders  were  allowed  to  enter  the  optical 
path,  the  experimental  results  would  be  negatively  affected.  No  calculation  was  performed 
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to  identify  the  diffraction  limited  spot  size  present  in  the  focal  plane  where  the  iris  is 
positioned;  the  iris  was  adjusted  visually. 
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Figure  3.12:  Exploded  view  of  lens  systems  on  Path  1  and  Path  2  of  the  experiment  design 
highlighting  the  translated  SLM  planes  and  lens  focal  lengths. 


The  second  purpose  of  this  4-/  lens  system  is  to  translate  the  SLM  plane.  This  plane 
is  located  at  the  focal  plane  of  the  second  lens  in  the  system.  This  plane  translation 
is  necessary  to  prevent  the  addition  of  phase  curvature.  The  4mm  Gaussian  amplitude 
filter  (GAF)  is  placed  at  this  point. 

Following  the  GAF  on  Path  i  is  a  second  4-/  lens  system  (LS2)  shown  in  Fig.  3.12 
which  serves  the  purpose  of  translating  the  SLM  plane  again  for  the  placement  of  the  LC 
variable  retarder  (VR).  The  lenses  used  are  plano-convex  and  have  a  100mm  focal  length. 
Prior  to  the  VR,  there  are  several  other  optical  elements  requiring  discussion.  A  HWP 
is  placed  within  the  lens  system  to  transform  the  vertically  polarized  beam  to  horizontal 
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polarization.  After  the  now  horizontally  polarized  beam  leaves  the  second  lens  in  the 
system,  it  passes  through  another  PBS  for  recombination. 

Following  the  SLM  on  Path  2  is  a  4-/  lens  system  (LS3)  shown  in  Fig.  3.12  which 
again  serves  the  same  purposes  as  the  first  4-/  system  on  Path  1.  This  lens  system  also 
has  an  iris  at  the  focal  plane  between  the  lenses.  The  lenses  used  are  plano-convex  with  a 
250mm  focal  length.  The  iris  again  is  used  to  remove  unwanted  diffraction  orders.  At  the 
focal  plane  of  the  second  lens  in  the  system  (the  translated  SLM  plane)  is  a  3mm  GAF. 

Following  the  GAF  on  Path  2  is  a  second  4-/  lens  system  (LS4)  shown  in  Fig.  3.12 
which  translates  the  SLM  plane  again  for  placement  of  the  previously  mentioned  VR.  The 
lenses  used  are  plano-convex  with  a  225mm  focal  length.  After  the  vertically  polarized 
beam  passes  through  the  second  lens  in  the  system,  it  reflects  through  the  PBS  previously 
mentioned  for  recombination. 

The  PBS  at  the  end  of  path  recombines  the  Path  1  and  Path  2  beam  components.  The 
final  element  the  recombined  beam  must  pass  through  is  the  VR,  placed  in  the  translated 
SLM  plane.  The  ThorLabs  LCC1223-A  full- wave  LC  VR  is  electrically  controlled  by  a 
ThorLabs  LCC25  LC  controller.  In  order  to  only  control  the  phase  of  the  incident  EGSM 
beam,  the  linearly  polarized  input  beam  must  have  a  polarization  axis  aligned  with  the 
optical  axis  of  the  VR  [22].  As  the  voltage,  Vrms,  is  increased  on  the  controller,  the  phase 
offset  in  the  beam  is  decreased.  The  VR  is  oriented  to  retard  the  horizontally  polarized  axis, 
i.e.,  alter  the  optical  path  length  of  Path  1  with  respect  to  Path  2  and  provide  ZBjcy.  Note 
the  VR  was  not  calibrated  due  to  time  constraints  and  was  set  to  minimize  the  amount  of 
retardance  applied  for  purposes  of  this  experiment. 

The  beam  exiting  the  VR  constitutes  an  instantaneous  realization  of  an  EGSM  beam 
at  the  source  plane.  Immediately  after  the  VR  is  a  1000mm  focal  length  plano-convex  lens. 
The  purpose  of  this  lens  is  to  focus  the  beam  to  the  sensor.  The  reason  this  lens  is  required 
is  due  to  fact  that  the  sensor  is  smaller  than  the  SEM.  Given  the  constraints  on  the  table 
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size,  there  was  not  room  to  feasibly  inelude  an  additional  4-/  lens  systems  to  demagnify  the 
beam  prior  to  the  sensor.  For  this  reason,  the  experiment  measures  the  EGSM  observation 
plane  rather  than  the  souree  plane. 

3.4.2  Polarization  Analyzer. 

To  eolleet  the  desired  data  from  the  generated  EGSM  souree,  the  beam  was  passed 
through  the  previously  mentioned  foeusing  lens  to  a  quarter- wave  plate  and  linear  polarizer. 
These  two  elements  eonstitute  what  is  referred  to  as  a  polarization  state  analyzer.  Following 
this  analyzer  is  the  deteetor,  plaeed  at  the  foeal  plane  of  the  1000mm  foeusing  lens.  This 
foeal  plane  is  the  EGSM  observation  plane. 
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Figure  3.13:  Sehematie  illustration  of  polarization  state  analyzer  eomposed  of  a  foeusing 
lens,  quarter-wave  plate,  linear  polarizer,  and  imaging  eamera. 


The  deteetor  used  is  a  Edmund  Opties  5012M  CMOS  imaging  sensor  with  a  2560  x 
1920  pixel  resolution  and  piteh  of  2.2yum  [5].  The  region  of  interest  for  purposes  of  this 
experiment  was  eropped  to  1024x  1024  to  reduee  the  amount  of  data  stored  in  the  eolleeted 
files.  The  images  were  able  to  be  eropped  due  to  the  faet  that  light  ineident  on  the  sensor 
was  not  seattered  outside  the  ehosen  region.  Additionally,  the  deteetor  was  not  ealibrated, 
whieh  eould  lead  to  possible  errors  in  the  gathered  data. 

3.4.3  Data  Collected. 

With  this  experimental  design,  the  desired  measurements  at  the  EGSM  observation 
plane  are  the  Stokes  parameters,  DoP,  and  SDoC.  None  of  these  measurements  are 
direetly  available  from  the  imaging  sensor  so  these  results  must  be  eomputed  from  the 
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irradiance  incident  on  the  sensor.  With  the  help  of  the  previously  mentioned  polarization 
state  analyzer,  the  QWP  and  LP  are  set  to  specific  orientations  as  detailed  in  Tab.  3.1.  The 
irradiance  images  gathered  for  each  set  of  orientations  allows  the  unnormalized  Stokes 
parameters  to  be  calculated  as  [20] 

S  0  =Ih  +  Iv 

S 1  =Ih  -  h 

(3.21) 

S  2  =1+45  -  1-45 

S3  =Ir-Il 

where  5  o  is  the  incident  irradiance,  S  i  is  the  horizontally  polarized  irradiance  In  minus 
the  vertically  polarized  irradiance  /y,  5  2  is  the  45-deg  polarized  irradiance  7+45  minus  the 
135-deg  polarized  irradiance  1-a5,  and  S3  is  the  right-hand  circularly  polarized  irradiance 
1r  minus  the  left-hand  circularly  polarized  irradiance  Ir-  The  wave  plate  and  polarizer 
orientations  required  to  collect  these  irradiance  components  are  defined  in  Tab.  3.1. 


Table  3.1:  Polarization  Analyzer  Orientations 


Irradiance 

QWP 

LP 

Ih 

0-deg 

0-deg 

Iv 

90-deg 

90-deg 

7+45 

45-deg 

45-deg 

1-45 

-45-deg 

-45-deg 

Ir 

-45-deg 

0-deg 

II 

45-deg 

0-deg 
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With  the  Stokes  parameters  defined  and  calculated,  the  DoP  for  the  total  beam  is  then 


calculated  as 


DoP  = 


2  +  52  +  52 


(3.22) 


and  the  degree  of  linear  polarization  (DoLP)  is  calculated  as 


DoLP  = 


’2  +  52 


(3.23) 


which,  if  5  3  «  0,  as  is  the  case  in  some  of  the  experimental  results,  DoP  «  DoLP. 

The  last  desired  measurement  is  the  SDoC.  This  is  not  directly  measurable  nor  able  to 
be  calculated  using  the  gathered  irradiance  images  because  an  electric  field  would  have  to 
be  captured  by  the  detector.  Given  the  degree  of  coherence  (DoC)  of  each  CSDM  element 
as  defined  by 


^apiPl^  P2) 


a  =  x,y 
--  x,y 


(3.24) 


V^ays(Pl  ’  Pl)^apiP2,  Pi) 
a  comparable  measurement  needs  to  be  taken  to  obtain  this  quantity.  Thus,  taking  the 

square  of  the  modulus  yields 

W,^(PUP2)W:^(PUP2) 


IPa/sl 


a  =  x,y 
--  x,y 


(3.25) 


5  a(Pl)S/3(p2) 

Using  the  gathered  irradiance  images  as  compared  to  the  electric  fields  and  applying  the 
Gaussian  Moment  Theorem  yields 


(Ia(Pl)IdP2))  ={E:(P,)E;(P2)EJp,)E^(P2)) 

={EI(p:)E^(P,)){E;(p2)E^(P2))  +  {El{p,)Ep{p2)){E;{p2)E^{p,)) 
=SMSp{P2)  +  ^apiPuPlWlpiPuPl) 

(a  =  x,y  and  /3  =  x,y) 


This  simplifies  to 


2  _  (■^ff(Pl)-^ys(P2)) 

S  aipl)S  pipl) 


(3.27) 
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which  is  the  normalized  fourth-order  eorrelation  funetion  (FOCF)  expanded  in  terms  of  the 
DoC  [4,  23], 


For  the  purposes  of  providing  a  eleaner  result  when  eomparing  the  gathered 
experimental  data  to  the  simulated  and  theoretieal  data,  this  relationship  is  then  rearranged 
as 


{laiPl)W2))  =  W,p{puP2)W:p{puP2)  +  SMSpiP2) 


(3.28) 


where  the  irradianee  eorrelation  funetion  on  the  left  side  of  the  equation  is  readily  available 
from  the  experimentally  gathered  irradianee  images  and  the  values  in  the  sum  on  the  right 
side  of  the  equation  are  readily  available  from  the  simulated  and  theoretieal  data. 
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IV.  Analysis  and  Results 


Chapter  4  presents  the  results  of  two  different  EGSM  source  generation  experiments. 
The  first  (Experiment  I)  was  an  elliptically  partially  polarized  EGSM  source  with  a  fully- 
populated  CSDM.  The  second  (Experiment  II)  was  a  linearly,  partially  polarized  EGSM 
source  with  the  off-diagonal  elements  of  the  CSDM  equal  to  zero.  Table  4.1  details  the 
desired  and  actual  EGSM  source  parameters  used  in  Experiment  I  and  II.  Table  4.2  details 
the  required  phase  screen  values.  The  relations  between  the  desired  EGSM  source  and 
phase  screen  parameters  forms  a  system  of  coupled  nonlinear  equations  which  can  not  be 
analytically  inverted  [2].  To  determine  the  phase  screen  parameters  from  the  desired  source 
parameters,  constrained  nonlinear  optimization  was  used  to  find  the  optimal  parameters. 
Eor  simulation  and  theory,  512  points  per  side  and  a  spacing  of  15//m  were  used  to 
discretize  the  fields  along  Paths  1  and  2  in  Eig.  3.11.  These  numbers  were  chosen  to  match 
the  BNS  Model  P5 12-0635  SEM.  A  wavelength  of  d  =  632. 8nm  was  assumed.  The  results 
for  Experiment  I  and  II  are  detailed  in  following  sections. 
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Table  4.1:  EGSM  Source  Parameters 


Experiment  I 

Experiment  II 

Parameter 

Desired 

Actual 

Parameter 

Desired 

Actual 

A, 

1.3 

1.3 

Ax 

1.3 

1.3 

A, 

1 

1 

Ay 

1 

1 

^Bjcy 

0 

0 

^Bjcy 

0 

0 

o-x  (mm) 

2.8 

2.8 

O-X  (mm) 

2.8 

2.8 

(Ty  (mm) 

2.1 

2.1 

CTy  (mm) 

2.1 

2.1 

6xx  (mm) 

0.40406 

0.42643 

6 XX  (mm) 

0.40406 

0.40406 

6yy  (mm) 

0.30305 

0.30972 

6yy  (mm) 

0.30305 

0.30305 

6xy  (mm) 

0.44447 

0.41705 

6xy  (mm) 

0.44447 

0.44447 

\Bxy\ 

0.15 

0.14942 

\Bxy\ 

0 

2.5513e-6 

Table  4.2:  EGSM  Phase  Screen  Parameters 


Experiment  I 

Experiment  II 

Parameter 

Value 

Parameter 

Value 

Uax  (mm) 

2.4 

4.0.  (mm) 

2.9 

4.0.  (mm) 

1.4 

4.0.  (mm) 

1.7 

3.9143 

C^0. 

5.0552 

3.1416 

^0. 

6.3124 

E 

1 

E 

0.6225 
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4.1  Experiment  I  Results 

Experiment  I  was  completed  using  10,000  realizations  to  generate  the  EGSM  source. 
Eigure  4.1  shows  the  experimental  results  for  the  normalized  Stokes  parameters  compared 
to  the  results  of  10,000  simulations  and  theory.  The  images  are  organized  such  that 
the  theoretical,  simulation,  and  experimental  results  are  along  the  columns — theoretical 
results  are  Eigs.  4.1(a),  4.1(d),  4.1(g),  and  4.1(j);  simulation  results  are  Eigs.  4.1(b), 
4.1(e),  4.1(h),  and  4.1(k);  experimental  results  are  Eigs.  4.1(c),  4.1(f),  4.1(i),  and  4.1(1). 
Each  row  of  images  in  Eig.  4.1  is  a  Stokes  parameter — 5o  are  Eigs.  4.1(a),  4.1(b),  and 
4.1(c);  Si  are  Eigs.  4.1(d),  4.1(e),  and  4.1(f);  S2  are  Eigs.  4.1(g),  4.1(h),  and  4.1(i);  S3 
are  Eigs.  4.1(j),  4.1(k),  and  4.1(1).  Eigure  4.2  shows  slices  of  the  Stokes  parameters  for 
additional  visualization  of  the  results.  The  plots  are  organized  such  that  the  theoretical, 
simulation,  and  experimental  curves  overlay  each  other  for  each  Stokes  parameter — S  0  is 
Eig.  4.2(a),  S 1  is  Eig.  4.2(b),  S2  is  Eig.  4.2(c),  and  S3  is  Eig.  4.2(d). 

Overall  the  experimental  Stokes  parameters  in  Eig.  4.1  match  the  shape  of  the 
simulated  and  theoretical  parameters  well,  showing  good  agreement.  Noticeable  with 
the  experimental  results  for  Si  and  S2  in  Eigs.  4.1(f)  and  4.1(i)  is  an  apparent  shift 
of  the  absolute  value  of  the  maxima  at  the  center  of  the  beam.  This  is  due  to  spatial 
registration  issues.  Because  the  Stokes  parameters  (other  than  So)  involve  difference  when 
calculating  using  the  experimental  data,  any  misregistration  artifacts  are  amplified  [20]. 
The  misregistration  visible  in  these  parameters  is  a  result  of  the  rotating  polarizers  used  to 
capture  the  different  polarized  irradiance  components. 
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(x,  y,  z  =  1  m) 


X  (mm) 

(i) 


Figure  4. 1 :  Experiment  I  Stokes  parameter  results  compared  with  simulation  and  theory. 
The  rows  are  So,  S  i,  52,  and  5 3,  respectively,  while  the  columns  are  the  theory,  simulation, 
and  experimental  results,  respectively. 
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Figure  4.2:  Experiment  I  Stokes  parameter  results  eompared  with  simulation  and  theory. 
The  theory,  simulation,  and  experiment  sliees  plotted  together  for  eaeh  of  (a)  So,  (b)  Si, 
(e)  S2,  and  (d)  S3. 
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Further,  the  simulation  and  theoretieal  results  visible  in  Fig.  4.2  do  not  mateh  as  well 
as  published  results  [2].  This  is  due  to  the  faet  that  the  SLM  is  eropping  the  beam  prior 
to  passing  through  the  GAF.  This  limitation  was  unavoidable  in  the  experimental  design 
beeause  there  was  not  enough  spaee  available  on  the  optieal  beneh  to  allow  for  the  beam  to 
be  magnified  to  pass  through  the  GAFs,  then  demagnified  to  eontinue  down  eaeh  path. 


Figure  4.3:  Experiment  I  degree  of  polarization  results  for  theory,  simulation,  and 
experiment. 


Figure  4.3  shows  the  DoP  of  the  on-axis  field  with  theoretieal,  simulation,  and 
experimental  results  plotted  together.  The  figure  shows  that,  in  general,  the  DoP  of  the 
field  on  axis  ehanges  on  propagation  and  the  experimental  and  theory  approaeh  unity.  The 
simulation  DoP  does  not  approaeh  unity  due  to  the  previously  mentioned  issue  with  the 
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GAF  being  too  large  for  the  SLM.  The  rate  at  which  the  experimental  DoP  changes  does 
not  match  well  with  simulation  or  theory.  This  is  most  likely  due  to  the  experimental  phase 
retardance  being  unknown  yet  still  generating  the  CSDM  cross-terms  with  a  populated 
value  for  \Bxy\. 

Figure  4.4  shows  the  experimental  results  for  the  irradiance  correlation  function 
compared  to  the  results  of  simulation  and  theory.  The  images  are  organized  such  that 
the  theoretical,  simulation,  and  experimental  results  are  along  the  columns — theoretical 
results  are  Figs.  4.4(a),  4.4(d),  4.4(g),  and  4.4(j);  simulation  results  are  Figs.  4.4(b),  4.4(e), 
4.4(h),  and  4.4(k);  experimental  results  are  Figs.  4.4(c),  4.4(f),  4.4(i),  and  4.4(1).  Each 
row  of  images  in  Fig.  4.4  is  a  different  irradiance  correlation  result — {Ixixi,y\)Ixix2,y2)) 
are  Figs.  4.4(a),  4.4(b),  and  4.4(c);  {Ix{xi,yi)Iy{x2,y2))  are  Figs.  4.4(d),  4.4(e),  and 
4.4(f);  {Iy{xi,yi)Ix{x2,y2))  are  Figs.  4.4(g),  4.4(h),  and  4.4(i);  {Iy{xi,yi)Iy{x2,y2))  are 
Figs.  4.4(j),  4.4(k),  and  4.4(1).  Figure  4.5  shows  slices  of  the  irradiance  correlation 
functions  for  additional  visualization  of  the  results.  The  plots  are  organized  such  that  the 
theoretical,  simulation,  and  experimental  curves  overlay  each  other  for  each  irradiance 
correlation  result— (4(.ri,yi)4(.r2,y2))  is  Fig.  4.5(a),  {Ix{xi,yi)Iy{x2,y2))  is  Fig.  4.5(b), 
{Iy{xuy\)lx{x2,y2))  is  Fig.  4.5(c),  and  {Iy{xi,yi)ly{x2,y2))  is  Fig.  4.5(d). 

The  results  in  Figs.  4.4  and  4.5  showing  the  irradiance  correlation  functions 
{Iaix\,yi)Ip{x2,y2))  show  very  good  agreement  between  the  experimental,  simulated,  and 
theoretical  data,  thus  validating  the  ability  to  control  the  coherence  properties  of  the  EGSM 
beam. 
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Figure  4.4:  Experiment  I  irradianee  eorrelation  funetion  results  eompared  with  simulation 
and  theory.  The  rows  are  <4(xi,yi)4(x2,y2)),  <4(jci,yi)4fe,y2)),  {Ix{xi,yi)h{x2,y2)), 
and  {Ix{x\,y\)Ixix2,y2)),  respeetively,  while  the  eolumns  are  the  theory,  simulation,  and 
experimental  results,  respeetively. 
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Figure  4.5:  Experiment  I  irradiance  correlation  function  results  compared  with  sim¬ 
ulation  and  theory.  The  theory,  simulation,  and  experiment  slices  plotted  together 
for  each  of  (a)  <4(xi,yi)4(x2,y2)),  (b)  {Ix{xi,yx)ly{x2,y2)),  (c)  {Iy{xi,yi)Ix{x2,y2)),  and 
(d)  {Iy{Xuy\)Iy{X2,y2))- 


47 


4.2  Experiment  II  Results 

Experiment  II  was  completed  using  1,000  realizations  to  generate  the  EGSM  source. 
Eigure  4.6  shows  the  experimental  results  for  the  Stokes  parameters  compared  to  the  results 
of  1,000  simulations  and  theory.  The  images  are  organized  such  that  the  theoretical, 
simulation,  and  experimental  results  are  along  the  columns — theoretical  results  are 
Eigs.  4.6(a),  4.6(d),  4.6(g),  and  4.6(j);  simulation  results  are  Eigs.  4.6(b),  4.6(e),  4.6(h),  and 
4.6(k);  experimental  results  are  Eigs.  4.6(c),  4.6(f),  4.6(i),  and  4.6(1).  Each  row  of  images  in 
Eig.  4.6  is  a  Stokes  parameter — So  are  Eigs.  4.6(a),  4.6(b),  and  4.6(c);  S  i  are  Eigs.  4.6(d), 
4.6(e),  and  4.6(f);  S2  are  Eigs.  4.6(g),  4.6(h),  and  4.6(i);  S3  are  Eigs.  4.6(j),  4.6(k),  and 
4.6(1).  Eigure  4.7  shows  slices  of  the  Stokes  parameters  for  additional  visualization  of  the 
results.  The  plots  are  organized  such  that  the  theoretical,  simulation,  and  experimental 
curves  overlay  each  other  for  each  Stokes  parameter — So  is  Eig.  4.7(a),  S 1  is  Eig.  4.7(b), 
S2  is  Eig.  4.7(c),  and  S3  is  Eig.  4.7(d). 

Overall  the  experimental  Stokes  parameters  in  Eig.  4.6  match  the  shape  of  the 
simulated  and  theoretical  parameters  well,  showing  good  agreement.  Noticeable  for  the 
experimental  results  for  S 1  in  Eig.  4.6(f)  is  the  same  apparent  shift  of  the  absolute  value 
of  the  maxima  at  the  center  of  the  beam  as  in  Experiment  I.  This  is  again  due  to  spatial 
registration  issues  [20].  The  results  for  S2  and  S3  are  negligible  because  the  source 
parameters  used  have  CSDM  terms  Wxy  =  Wyx  =  0.  Eurther,  the  simulation  and  theoretical 
results  visible  in  Eig.  4.7  again  do  not  match  as  well  as  published  results  [2].  This  is  due  to 
the  fact  that  the  SEM  is  cropping  the  beam  prior  to  passing  through  the  GAE  as  described 
earlier. 
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Figure  4.6:  Experiment  II  Stokes  parameter  results  compared  with  simulation  and  theory. 
The  rows  are  5o,  ^  i,  52,  and  5 3,  respectively,  while  the  columns  are  the  theory,  simulation, 
and  experimental  results,  respectively. 
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Figure  4.7:  Experiment  II  Stokes  parameter  results  compared  with  simulation  and  theory. 
The  theory,  simulation,  and  experiment  slices  plotted  together  for  each  of  (a)  So,  (b)  Si, 
(c)  S2,  and  (d)  S3. 
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Figure  4.8  shows  the  DoP  of  the  on-axis  field  with  theoretieal,  simulation,  and 
experimental  results  plotted  together.  The  figure  shows  that,  in  general,  the  DoP  of  the 
field  on  axis  ehanges  on  propagation  and  the  experimental  and  theory  approaeh  unity.  The 
simulation  DoP  does  not  approaeh  unity  again  due  to  an  issue  with  the  GAF  being  too  large 
for  the  SLM.  The  rate  at  whieh  the  DoP  ehanges  shows  very  good  agreement  between 
experimental,  simulation,  and  theoretieal  results,  thus  validating  the  ability  to  eontrol  the 
polarization  properties  of  the  EGSM  beam. 


Figure  4.8:  Experiment  II  degree  of  polarization  results  for  theory,  simulation,  and 
experiment. 


Figure  4.9  shows  the  experimental  results  for  the  irradianee  eorrelation  funetion 
eompared  to  the  results  of  simulation  and  theory.  The  images  are  organized  sueh  that 
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the  theoretical,  simulation,  and  experimental  results  are  along  the  columns — theoretical 
results  are  Figs.  4.9(a),  4.9(d),  4.9(g),  and  4.9(j);  simulation  results  are  Figs.  4.9(b),  4.9(e), 
4.9(h),  and  4.9(k);  experimental  results  are  Figs.  4.9(c),  4.9(f),  4.9(i),  and  4.9(1).  Each 
row  of  images  in  Fig.  4.9  is  a  different  irradiance  correlation  result — {hix\,yi)Ixi^2^y2)) 
are  Figs.  4.9(a),  4.9(b),  and  4.9(c);  {Ix{xi,yi)ly{x2,y2))  are  Figs.  4.9(d),  4.9(e),  and 
4.9(f);  {Iyixuy\)lx{x2,y2))  are  Figs.  4.9(g),  4.9(h),  and  4.9(i);  {Iyixuyi)ly{x2,y2))  are 
Figs.  4.9(j),  4.9(k),  and  4.9(1).  Figure  4.10  shows  slices  of  the  irradiance  correlation 
function  for  additional  visualization  of  the  results.  The  plots  are  organized  such  that  the 
theoretical,  simulation,  and  experimental  curves  overlay  each  other  for  each  irradiance 
correlation  result— {I x(xi,yi)Ix(x2,y2))  is  Fig.  4.10(a),  {Ix(xuyi)Iy(x2,y2))  is  Fig.  4.10(b), 
{Iy(xi,yi)Ix(x2,y2))  is  Fig.  4.10(c),  and  {Iy(xuyi)Iy(x2,y2))  is  Fig.  4.10(d). 

The  results  in  Figs.  4.9  and  4.10  showing  the  irradiance  correlation  functions 
{Ia{x\,y\)Ij}ix2,y2))  show  very  good  agreement  between  the  experimental,  simulated,  and 
theoretical  data,  thus  further  validating  the  ability  to  control  the  coherence  properties  of  the 
EGSM  beam. 
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Figure  4.9:  Experiment  II  irradiance  correlation  function  results  compared  with  simulation 
and  theory.  The  rows  are  <4(xi,yi)4(x2,y2)),  <4(jci,yi)4fe,y2)),  {Ix{xi,yi)h{x2,y2)), 
and  {Ix{x\,y\)Ixix2,y2)),  respectively,  while  the  columns  are  the  theory,  simulation,  and 
experimental  results,  respectively. 
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Figure  4.10:  Experiment  II  irradiance  correlation  function  results  compared  with 
simulation  and  theory.  The  theory,  simulation,  and  experiment  slices  plotted  together 
for  each  of  (a)  <4(xi,yi)4(x2,y2)),  (b)  </x(-^i,yi)/>.(x2,y2)),  (c)  {Iy{xi,yi)h{x2,y2)),  and 
(d)  {Iy{Xuy\)Iy{X2,y2))- 
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V.  Conclusions  and  Recommendations 


5.1  Conclusions  of  Research 

The  experimental  design  demonstrated  in  this  research  has  been  shown  to  successfully 
produce  EGSM  beams.  The  generation  of  the  beam  was  accomplished  using  two  mutually 
orthogonal  polarized  components  of  initially  coherent  beams  with  Gaussian  amplitudes 
which  are  phase  modulated  by  a  sequence  of  two  mutually  correlated  random  phase 
screens  at  the  source  plane.  Coherence  and  polarization  statistics  were  calculated  on  the 
integrated  results  and  compared  to  published  simulated  and  theoretical  results  to  validate 
the  experimental  design  [2,  17]. 

Although  EGSM  beams  were  successfully  created,  the  accuracy  of  the  results  suffers 
without  the  calibration  of  all  electro-optical  components  used  in  the  beam  generation  and 
data  capture.  However,  based  on  the  results  obtained  in  this  research,  the  experimental 
design  provides  a  viable  method  for  generating  EGSM  beams  with  prescribed  coherence 
and  polarization  properties.  With  the  use  of  EC  SEMs,  exchangeable  amplitude  filters, 
and  a  EC  VR,  the  design  is  highly  configurable  to  generate  EGSM  beams  with  a  wide 
range  of  properties  for  application  in  free-space  optical  communication,  remote  sensing, 
and  directed  energy. 

5.2  Recommendations  for  Future  Research 

In  this  research,  the  experimental  design  was  created  and  tested  for  viability.  Euture 
work  would  first  entail  a  more  thorough  effort  on  calibration  of  the  electro-optical  devices 
used,  i.e.,  the  SEMs,  the  VR,  and  the  detectors.  This  would  both  reduce  error  associated 
with  aberrations  and  linear  phase  response,  but  also  allow  the  ability  to  accurately  measure 
all  spectral  Stokes  parameters  to  provide  more  accurate  results. 
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In  addition  to  calibration,  new  SLMs  and  detectors  could  be  used.  The  new  SLMs 
could  provide  a  higher  resolution,  greater  pixel  density,  larger  array  size,  16-bit  operation, 
and  compatibility  with  64-bit  host  operating  systems  to  improve  memory  shortcomings. 
A  division  of  amplitude  polarimeter  could  be  used  as  a  new  detector,  allowing  the 
simultaneous  capture  of  three  full  resolution  polarimetrically  filtered  images  from  which 
Stokes  parameters  can  be  computed  [20].  This  would  eliminate  the  need  for  a  rotatable 
polarization  analyzer  which  may  cause  image  registration  errors. 

With  regard  to  the  SLM,  this  device  could  be  used  for  more  than  simply  displaying  the 
phase  screens.  There  is  potential  for  applying  the  amplitude  filter  and  phase  retardance  on 
the  SLM.  Combining  these  elements  on  the  SLM  would  eliminate  the  need  for  the  GAFs, 
VR,  and  two  of  the  four  4-/  systems  in  the  experimental  design.  The  removal  of  these 
elements  would  reduce  error  associated  with  aberrations  or  misalignment.  It  would  also 
provide  a  greater  potential  of  producing  non-Gaussian  electromagnetic  Schell-model  (SM) 
beams  or  “designer  beams.” 

Lastly,  with  regard  to  the  physical  design  of  the  experiment,  the  4-/  lens  systems  could 
be  configured  to  demagnify  the  SLM  plane  so  the  EGSM  source  plane  could  be  captured  by 
a  detector  rather  than  the  EGSM  observation  plane.  Additionally,  the  EGSM  source  plane 
could  be  simulated  to  propagate  through  the  atmosphere  using  either  an  atmospheric  phase 
wheel  or  using  additional  SEMs  to  display  atmospheric  phase  screens.  Eurther,  the  source 
plane  could  be  reflected  off  an  object  to  analyze  the  scattered  light.  Both  of  these  proposals 
have  a  direct  impact  on  experimentally  verifying  theory  using  EGSM  beams  for  free-space 
optical  communication  and  remote  sensing. 
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